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We consider the diffusion and spreading of chainlike molecules on solid 
surfaces. We first show that the steep spherical cap shape density profiles, ob- 
served in some submonolayer experiments on spreading polymer films, imply 
that the collective diffusion coefficient Dc(9) must be an increasing function 
of the surface coverage 6 for small and intermediate coverages. Through simu- 
lations of a discrete model of interacting chainlike molecules, we demonstrate 
that this is caused by an entropy-induced repulsive interaction. Excellent 
agreement is found between experimental and numerically obtained density 
profiles in this case, demonstrating that steep submonolayer film edges nat- 
urally arise due to the diffusive properties of chainlike molecules. When the 
entropic repulsion dominates over interchain attractions, Dc{0) first increases 
as a function of 6 but then eventually approaches zero for — > 1. The max- 
imum value of Dc{0) decreases for increasing attractive interactions, leading 
to density profiles that are in between spherical cap and Gaussian shapes. We 
also develop an analytic mean field approach to explain the diffusive behav- 
ior of chainlike molecules. The thermodynamic factor in Dc(9) is evaluated 
using effective free energy arguments, and the chain mobility is calculated 
numerically using the recently developed dynamic mean field theory. Good 
agreement is obtained between theory and simulations. 

PACS numbers: 68.35.Fx, 68.10.Gw, 83.10.Nn 



I. INTRODUCTION 

The dynamics of polymers on solid surfaces is an interesting theoretical problem that 
has important applications related to thin surface films. A central role in such systems is 
played by the diffusive dynamics of the chains. While the diffusion of adatoms and small 
molecules on surfaces is an extensively studied problem , there are relatively few studies 
of the diffusion of more complicated molecules, and especially polymers on surfaces 
The pioneering experiments of George et al. || on short-chain alcanes on metal surfaces 
revealed that the coverage dependence of the collective diffusion coefficient Dc{6) may show 
unusual features. For some molecules, it is approximately constant while for others, it is an 
increasing function of coverage for a wide range of values of 9. 

There has also been substantial interest recently on the spreading dynamics of molecu- 
larly thin oil films on solid substrates [|6|-4T(i|j. Many of these experiments are performed by 



depositing tiny, very flat droplets on surfaces. In the limit where the film becomes less than 
one monolayer in thickness, the spreading molecules are all in contact with surface forming 
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a 2D molecular gas. In this regime, the dynamics of spreading and consequently the density 
profiles of the film edges are determined by the diffusion of the molecules. Experimental 
studies of spreading in the submonolayer regime reveal that most of the measured film 
edge profiles are not well approximated by the Gaussian function but assume a steeper shape 
fTT| that can be well fitted by a spherical cap in the rotationally symmetric case, i.e. droplet 
spreading f|||. 

In the regime where the diffusion equation (Fick's law) is a valid description, a non- 
Gaussian profile indicates that the collective diffusion constant Dq{9) must have nontrivial 
coverage dependence. In the present context this was first realized in conjunction with the 
experiments of Albrecht et al. [Q, where the observed spherical cap shaped droplets could 
be reproduced by Monte Carlo simulations of spreading of flexible chains in 2D. Analyzed 
in terms of Fick's law with a coverage-dependent D c (6), it was concluded in Ref. [|11[] that 
D c {9) is a strongly increasing function of coverage up to 9 « 0.4 because of entropic repulsion 
between the chains. 

In this work, our aim is to extend the work presented in Ref. [|12|] to carry out a system- 
atic study of the spreading and diffusion of interacting chainlike molecules on surfaces. To 
begin with, we present an analysis of density profiles of spreading droplets with qualitatively 
different forms of D c {9). Using Fick's law, we show how the spherical cap type of shapes 
result from the increase of D c (6) vs. 9, and then quantitatively determine D c {9) from the 
experimental profiles of Ref. || up to 9 « 0.4. By extensive Monte Carlo simulations of a 
discrete model of interacting chainlike molecules, we then compute D c (8) for all coverages, 
and demonstrate that the chainlike nature of the molecules induces a strong repulsive inter- 
action. When this dominates over interchain attractions, Dc{9) first increases as a function 
of 9 as deduced from the experiments, but eventually must go to zero for 6 — > 1. The rela- 
tive maximum of Dc{9) decreases for increasing attractions. We also examine the profiles of 
spreading droplets in detail for various interactions. For the case of pure entropic repulsion, 
we find excellent agreement between experimental and numerically obtained profiles. This 
demonstrates that steep density profiles can be obtained from energetic considerations with- 
out having to assume that the film edge acts as a phase boundary between a 2D condensate 
and a vapor phase, as suggested previously ||. For increasing attractive interactions, the 
corresponding droplet shapes are in between spherical cap and Gaussian shapes. 

Finally, we develop an analytic mean field approach to explain the behavior of the collec- 
tive diffusion coefficient for chainlike molecules by using the fact that Dc{9) can be expressed 
as the product of thermodynamic factor Xo 1 and the mobility m. The thermodynamic fac- 
tor is evaluated using effective free energy arguments, and the chain mobility is calculated 
numerically using the recently developed dynamic mean field theory [O. Good agreement 
is found between theory and the numerical simulations. 

The outline of the paper is as follows. In Sec. II we describe the analysis of spreading 
profiles from the diffusion equation which gives partial information about the collective 
diffusion coefficient. The lattice model used in this work is explained in Sec. III. The 
results of numerical simulations are presented in Sec. IV. We discuss briefly the results for 
athermal chains presented in Ref. ||12|| , and present complete results for chains with attractive 
interchain interactions. In Sec. V we present the mean field theory for the thermodynamic 
factor Xo 1 an d the mobility m. In Sec. VI we use the chain-chain pair distribution function 
to extract the effective pair interaction potential for athermal chains. Finally, in Sec. VII 
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we briefly summarize our main results. 



II. ANALYSIS OF SPREADING PROFILES FROM THE DIFFUSION 

EQUATION 

As discussed in the Introduction, experimental measurements of spreading density pro- 
files of a polymer film in the submonolayer regime can be used to obtain information about 
the coverage dependence of the collective diffusion coefficient Dc{9). This was first realized 



by Herminghaus et al. [ I I ] who considered the spreading of polydimethylsiloxane (PDMS) 
on metal surfaces, including the data of Ref. ||. A typical late-time spreading profile taken 
from this experiment is shown in Fig. 1. The height of the profile is measured in Angstroms 
(A), but it actually describes the average areal density distribution of polymers on the sur- 
face. As noted in Ref. |§, this and the other submonolayer profiles can be fitted very well 
with spherical cap function 8{x) = 9q + \/r 2 — x 2 which is shown in Fig. 1 with a dotted line 
(#o and r are fitting parameters). In the same figure we also show a Gaussian profile with 
a dashed line for comparison. While the shape of the profile clearly follows the spherical 
cap shape for most of the range of densities, at the lowest coverages the film edge displays 
a Gaussian shaped tail. 

A qualitative way of analyzing the implications of the density profiles for Dc(0) is to 
apply the non-linear diffusion equation Jl],|lT| 



dO(x,t) _ d_ 
dt dx 



Dc{d) ^x~ 



(1) 



Starting from a delta function type of initial profile, the solution 9(x,t) of Eq. (1) is 
an exact Gaussian profile for all times when Dq = const. In Fig. 2 we show typical results 
from the numerical solution of Eq. (1) for three choices of Dc{0), including a monotonously 
decreasing and increasing functions of coverage Dc(0) = 1—0.990 and a Dc(0) = 0.01+0. 99$, 
respectively. For the first case, profiles are obtained that are below the Gaussians for high 
and intermediate coverages. On the other hand, for cases where Dc{0) increases, profiles 
are above the Gaussian solution. This is exactly what was observed in the experiments (see 
Fig. 1) and thus we can conclude that Dc(9) must be an increasing function of coverage for 
the range of covereges corresponding to the density profiles (see also Ref. [|TT|j). 

It is possible to determine Dq{9) more quantitatively by solving density profiles from Eq. 
(1) and matching them with the experimental ones. We did this by using a fitting function 
for D c {9) of the form D c {9) — c\ + c 2 tanh[c 3 (6' — c 4 )], where c±, c 2 , c 3 , and c 4 are fitting 
parameters. By adjusting these parameters it is relatively easy to obtain profiles with shapes 
closely matching the experimental ones of Ref. ||, with c\ = 0.159, c 2 = 0.144, C3 = 5.174, 
and C4 = 0.345. The resulting profile is shown in Fig. 1 with a solid line. The matching 
of the profiles also fixes the normalized coverage scale from Fick's law corresponding to the 
effective layer thickness in the experiments. The analysis reveals that Dc{6) at very low 
concentrations is almost constant and then monotonically increases up to about 9 = 0.4, in 



agreement with Ref. [11|. The resulting Dq{9) from the best fits is shown in the inset of 
Fig. 1. We not that the behavior of Dq{9) at higher values of 9 cannot obtained from this 
procedure, because the maximum value of the experimental profile at hand is already less 
than 9 = 0.4. 
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III. MODEL OF CHAINLIKE MOLECULES 



To determine Dc{0) for all values of 9 and to study the effects arising from the chainlike 
nature of the diffusing molecules on surfaces, we will in the present work use a discrete model 



called the fluctuating bond (FB) model [13]. This minimal model of polymers is widely used 
for simulations of many-chain systems [|TjJ. The idea in the FB model is that it is a coarse- 
grained model of real polymer chains. Real polymers, such as simple polyethylene consist of 
repeated segments of CH 2 monomers where carbon atoms are bound to each other forming a 
long chain. The bond length and the bond angle between adjacent carbon atoms are almost 
fixed. However, the torsional angle between adjacent bonds can have different values, and 
thus the end-to-end distance of a long chain of CH 2 segments may vary significantly. On a 
coarse-grained level then, such a chain can be described by a reduced number of effective 
segments. 

In the 2D version of FB model the chains in the model consist of connected segments that 
occupy sites on a square lattice. Each segment prohibits all other segments from occupying 
its nearest or next nearest neighbor lattice sites. In the model (see Fig. 3), the distance 
between adjacent segments i can vary between 2 < i < vT3 in lattice units, where the upper 
limit prevents bonds from crossing each other. With these restrictions there are 36 different 
bond vectors in the model. The FB model has been shown to give static properties, such as 



pressure, in full agreement with simulations of continuum models [16|. Furthermore, the FB 



model incorporates the same type of dynamics for a single chain as the continuum Rouse 



model 113 



To consider the general case where there are attractive interactions between the chains 
and the flexibility of individual chains can also vary, we have used the following Hamiltonian: 



k B T k B T 



\ 12 / 

a \ a 



rp n N FB — 1 

" E E (2) 



In Eq. (2), the first term is a Lennard- Jones type of potential, where J = Ae/ksT is 
the strength of the interaction, n,2';i,»' is the distance between segments i and i' of different 
chains I and and a was chosen such that the potential minimum was at the distance of 
two in lattice units. The cut-off radius of the potential was v'lO lattice units, after which 
the potential is small enough to give a negligible contribution. In the summations, n is 
the number of chains and Np B is the number of segments in each chain. The second term 
controls the stiffness of each chain, with <pi^ being the angle between two adjacent bonds i 
and i + 1 in chain I and K = E s /kBT the stiffness parameter. For increasing E s or decreasing 
temperature T, the chains become more stiff. 

Dynamics is introduced in the model by Metropolis moves of single segments, with a 
probability of acceptance min[e~ AWoff//A:sT , 1], where A7i e ff is the energy difference between 
final and initial configurations for acceptable moves to a nearest neighbor site, for which site 
exclusion and bond length restrictions must be satisfied. One Monte Carlo (MC) time step 
is defined as an attempt to move each segment of every chain. In Fig. 3 we show a typical 
configuration of a chain in the FB model for Np B — 6. One possible move of a segment is 
shown by the dashed line. It should be noted that since the dynamics within the FB model 
consists of single segment moves only and there are no direct translational modes, the limit 
of a rigid rod is not well defined p 
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IV. NUMERICAL SIMULATIONS 



A. Athermal Chains 

To quantitatively determine the coverage dependence of diffusion, we have performed 
extensive MC simulations using the FB model explained in Sec. III. First, we consider 



the case of fully flexible, athermal chains for which J = K = |fL2|| . The linear size of 
the 2D square lattice we used was L = 180 for most cases. To calculate D c {0) we used 
the temporal decay of the Fourier transformed density autocorrelation function S(k,t) = 
S{k,0)e~ k2Dc ^ t , where S(k, 0) becomes constant in the hydrodynamic limit k — > 0. The 
density-fluctuation autocorrelation function is defined as S(r,t) = (59(r,t)59(0,0)), where 
59(r, t) = 9(r, t) — (9(r, t)), and S(k, t) is its Fourier transform fll8| . It is calculated separately 
for each fixed value of the average normalized coverage 9 = (9(r, t)) which for the FB model 
is defined to be 9 = AuNfb/L 2 . In Fig. 4 we show results of these calculations for the case 
Nfb = 6 with circles. The results are normalized with Di, which is the diffusion constant 
of a single monomer (Npg — 1) in the limit 9 — > 0. Initially, Dq{9) increases up to 9 m 0.7, 
after which it rapidly approaches zero. The initial behavior of Dc{9) is in good agreement 
with results obtained from the shape of the spreading profiles. 

It is interesting to compare the behavior of Dc{9) with the tracer diffusion coefficient of 
individual chains Dt{9) defined by 

1 n 

Dr = }™^um)-m\ 2 ), (3) 

where fi(t) is the position vector of the 2 th chain at time t. In Fig. 5 we show the results 
from simulations for the athermal case with circles. In the limit 9 — > Dc(0) = Dt(0), and 
thus we use the same normalization as in Fig. 4. As expected from increased interchain 
blocking, Dt{9) is a monotonically decreasing function of 9. This is a striking example of 
how fundamentally different the behavior of the two types of diffusion coefficients can be. 

We have also studied the effect of chain length and stiffness to collective diffusion, and 



a summary of the results can be found in Ref. [12]. When the length of the chains increases 



diffusion slows down, but the relative maximum of Dc(9) becomes more pronounced. With 
increasing stiffness diffusion slows down, too, but now the maximum of Dq{9) becomes less 
pronounced. 



B. Chains with Attractive Interactions 

To study the effect of attractive interactions between polymers we used the FB model 
with the Hamiltonian of Eq. (2) with two interaction parameters, namely J = —0.5 and 
J = —1.0. For the results presented here, we considered fully flexible chains (K = 0) of 
length Nfb = 6. The sizes of the lattices in the MC simulations varied from 100 x 100 to 
180 x 180. 

The results for the two CclSGS clS obtained from the decay of S(k,t) are shown in Fig. 4. 
As the strength of the attractive interaction increases, diffusion slows down but the overall 
behavior as a function of 9 is qualitatively similar to the athermal case. The main influence 
of the attraction is to significantly reduce the relative height of the maximum of Dq{9) ||19|| . 
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To compare with the athermal case, we also calculated Dt{0) for which the results are 
shown in Fig. 5. When the strength of the interaction increases, Dt{9) decreases more 
rapidly function of concentration. 



C. Spreading of Droplets 

We also used the FB model to directly simulate the spreading of 2D submonolayer 
droplets. The chains were initially confined to a circular region, and after equilibration the 
spatial constraints were removed. The consequent spreading was monitored and correspond- 
ing density profiles calculated as a function of time. We find that within the accuracy of 
the fits, there is a linear relationship between the experimental and MC time scales which 
supports the use of the single- segment dynamics within the FB model. 

In Fig. 6 we show a comparison between three experimental submonolayer profiles |§, 
profiles obtained from the MC simulations for athermal chains, and profiles computed using 
Eq. (1) with the tanh fitting function for Dc(9) as explained previously. The agreement 
between all the profiles is excellent demonstrating the consistency of our approach. It also 
shows that the behavior of the PDMS polymers used in the experiment can be most simply 
explained in terms of athermal chain dynamics, with the entropic repulsion dominating. 
This is in contrast to the experiments of Ref . || , where the steep film edges were assumed 
to be a ID phase boundary between a 2D condensate and a vapor phase. Our results here 
demonstrate that no such assumptions need to be made; the steep spherical cap shapes of 
submonolayer droplets are expected to be a generic signature of strongly repulsive effective 
interactions. 

In Figs. 7(a) and 7(b) we show additional spreading profiles for the cases J = —0.5 
and J = —1.0, respectively. The changes seen in Dc{0) in Fig. 4 can also be seen in the 
corresponding spreading profiles. In Figs. 7(a) and 7(b) the simulation results are shown 
with circles, a spherical cap fit with a solid line and a Gaussian fit with a dotted line. The 
main result is that with increasing attractive interactions, the shape of profiles changes from 
the spherical cap towards a Gaussian shape (except for the highest coverages). These results 
show the intimate connection between Dc{0) and spreading profiles; in the regime where 
Dc(0) is almost constant, the corresponding profile shape is close to the Gaussian limit. 



V. MEAN FIELD THEORY FOR COLLECTIVE DIFFUSION 

To better understand the somewhat unusual behavior of collective diffusion of chainlike 
molecules, we start from the Green-Kubo relation Jl| 

D c = lim _— / dt(J(0) ■ J(t)) , (4) 
t^oo 2{(on) z ) Jo 

where n is the number of chains, J(t) = Z)£=i^i(0 is the total particle current, and 
{(Sn) 2 ) the mean square fluctuation of chains (in a finite area A). In terms of the mean 
square displacements of the individual chains, this can be written as 
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= Xo l m, (5) 

where r^(i) is the position of the center of mass of the i th chain at time t. In this 
equation, the term Xo 1 — ( n ) / (($ n ) 2 ) defines the thermodynamic factor, and the remainder 
of the equation defines the mobility m. The thermodynamic factor is related to the density 
fluctuations of the system, while the mobility can be written as 

^ = i lim 3 -((Af CM ) 2 ), (6) 

where Ar C M = ^cuit) — r*c M (0) is the displacement of the center of mass of all the 
chains, and r C M(f) = (V n ) 5Z"=i ^KO- I n the theory presented here, we will treat the two 
factors Xo 1 and m separately. 



A. Thermodynamic Factor 

To estimate the thermodynamic factor we consider a generalization of the simple ther- 



modynamic theory presented for athermal chains in Ref. [|T^]. We take as a starting point 
an effective Helmholz free energy F as 

F = F - na9 b - k B Tln[————] - nk B T\nw. (7) 

(M — ny.nl 

The first term F is a constant, while the second term E int = na6 h comes from attractive 
interactions between segments of different chains and is temperature dependent. We 
calculated E int directly from the MC simulations at different coverages to estimate a and b 
and verified that this approximation is well satisfied. For the cases J = —0.5 and J = —1.0 
the resulting values of these parameters are a/ksT = 1.8, b = 1.6, and a/ksT = 2.7, b = 1.4, 
respectively, while for the athermal case a = b = 0. The third term comes from the entropy 
of the center of the mass of the n polymers and it is here approximated by the expression 
for a 2D Langmuir gas on a lattice [Q. The parameter M = L 2 denotes the total number 
of lattice sites and thus 9 = nN/M (for the FB model 9 = AnNpg/M due to the exclusion 
rules). The last term in Eq. (7) is the entropic contribution from the chainlike nature of 
the molecules, where w is the number of possible configurations of each chain, and is a 
model dependent quantity. For the present case of chainlike molecules, we approximate it 
by decoupling the total number of configurations into the product of two terms 

w{9) = wl(9)w^- 2 (9), (8) 

where W\{9) is the entropy arising from a segment at the end of the chain, and W2{9) 
from each segment in the middle of the chain. For the FB model with different interactions, 
we have numerically determined w\(9) and W2{9). In Fig. 8(a) we show the behavior of w\ 
and u>2 with three different values of J, namely J = 0, J = —0.5 and J = —1.0. These 
quantities can be easily interpolated for all coverages [|(J . 



The chemical potential can be calculated from F by using /z = (dF/dn)xy which gives 
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where p = 8/N is the number of chain molecules per unit area with N segments, and 
a' = a/ksT. It can be shown that the thermodynamic factor can be written as Xq 1 = 
6[d(p/k B T)/d6] [§, and thus we obtain 

X0 - 1 = 1 - o ' K 6 + lK + ^- e ^M. (10) 

In Fig. 8(b) the markers show results for Xo 1 ^) as obtained from accurate MC sim- 
ulations of the density fluctuations of the FB model (with Nfb = 6) in equilibrium, as 
extrapolated to an infinite system size. At very high concentrations density fluctuations are 
so small that it is very difficult to obtain accurate results. Furthermore, with attractive 
interactions (J 7^ 0) for coverages 9 > 0.8 the dynamics of the system slows down |Tj| and 
thus we here present results only for smaller coverages. Using the approximation of Eq. (|10"D 
with numerically determined w\{9) and 102(6) of Fig. 8(a) with Nfb = N for the ather- 
mal case, the results show that the true magnitude of Xo 1 ^) is somewhat underestimated 
throughout the range of coverages However, if the magnitude of the thermodynamic 
factor is known for some coverages, the effective chain length N appearing in Eq. (|8]) can 
be used as an additional fitting parameter to improve the results. In Fig. 8(b) we show the 
results of this approach, with only one parameter N fitted to our MC data for Xo The 
corresponding values of the parameters for J = —0.5 and J = —1.0 are N = 11, and N = 6, 
respectively. For the athermal case N = 1Q gives the best results. Thus, for increasing 
attraction, the approximation of Eq. (H) seems to become more accurate. 



B. Calculation of Mobility 

While the thermodynamic factor contains information about the equilibrium density 
fluctuations, the mobility m is determined by the dynamics of the center-of-mass motion of 
the particles. To calculate m theoretically we use the recently developed Dynamical Mean 
Field (DMF) theory ]13|, which yields an approximate expression for m as 



£ 2 

m DMF = -jT, (11) 

where I is the effective jump length and T is the average jump rate. This formulation 
makes it very efficient to numerically evaluate m, and has recently been shown to give a 



very good approximation of the true m for various strongly interacting systems [T3| . 

In Fig. 9(a) we show mobilities calculated from MC simulations of m(8) for the FB model 
for the cases J = 0, J = —0.5 and J = — 1.0, using the definition of Eq. (6). The data are 
normalized with the mobility of one segment mi. In the same figure we show also the results 
calculated from the DMF theory of Eq. (|TTJ) . The effective jump length i has been estimated 
from the zero coverage limit, where Xo(0) = 1 and thus mi = D\. In the athermal case the 
DMF theory is in good agreement with simulation results, but with increasing attraction 
between the chains it starts to deviate more from the MC simulation results. This behavior 
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is physically reasonable, because attractive interactions strengthen the effect of dynamical 
correlations that are not included in the DMF theory |L3|. Despite this, DMF gives the 
qualitative behavior of m(9) rather well even in case of attractive chains. 

A commonly used method to approximate the mobility is based on the Darken equation 
|Ij]. It states that the mobility can be approximated by the tracer diffusion coefficient, i.e. 
m(9) ~ Dt(0). In Fig. 9(b) we show the complete results from simulations of Dt(0) with 
the interaction parameters J = 0, J = —0.5 and J = —1.0. In same figure, the markers 
denote the results of direct MC simulations of m(8). As can be seen from the comparison, 
for the chainlike molecules Dj-(9) and m(9) behave quite differently; for athermal chains 
even the curvatures of the two functions have opposite signs. 

VI. EFFECTIVE POTENTIAL 

The fact that the initial increase of Dc(0) is due to effective repulsive interactions between 
the molecules can be seen in Eq. (10), where the important term 1 — 9[dln(w)/d9] comes 
from entropic origin of the chainlike molecules [0,0]- With increasing attraction, this 
entropy- induced repulsion is compensated by the attractive E int , and the maximum value 
of D c {6) is reduced in magnitude. 

To study the effective potential corresponding to the entropic repulsion we calculated 
the chain-chain pair distribution function g(r) = J2u'(d( r — r i)d( r i')) HH; where is the 
position of the center of mass of chain i and prime indicates that terms with i = i' are 
to be omitted. From g(r) we extracted numerically the effective pair interaction potential 
V e (r)/kBT for athermal chains. This can be done by first calculating the direct correlation 
function c(r) from the Ornstein-Zernike relation |22| 



h(r) = c(r) + 9 J h(r')c(\r - r'\)dr', (12) 

where h(r) = g(r) — 1 and r is distance in lattice units. When c(r) is known, V^(r)/fceT 
can be calculated by using the hypernetted chain theory p3[ 



c(r) = -V e {r)/k B T + h{r) - log[/i(r) + 1]. (13) 

In Fig. 10 we show V e (r) jksT for the athermal case and for 9 = 0.25, which shows a 
strong repulsion extending up to several lattice cites |[24| . As a comparison in the same figure 
there is also a typical Lennard- Jones potential which is more repulsive at small distances. 
An interesting result of the analysis is that for athermal chains, all the computed pair 
correlation functions for coverages 9 < 0.7 and for several chain lengths collapse to a 
single function which is given by g{f) = G(r9 a /Np B ), where a w 0.38 and (3 w 0.55. These 
scaled correlation functions are shown in the inset of Fig. 10. 



VII. DISCUSSION AND CONCLUSIONS 

In this paper we have presented a systematic study of diffusion and spreading of chainlike 
molecules, in part inspired by the non-Gaussian submonolayer film profiles observed in most 
experiments. Using Monte Carlo simulations with the fluctuating bond model, we have 
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calculated the diffusion coefficients as a function of coverage, generalizing the results for 



athermal chains of Ref. |12j to chains with attractive interchain interactions. Typically, 
the collective diffusion coefficient Dq{9) increases initially and displays a maximum around 
9 ~ 0.7. The strength of the peak decreases with increasing attraction. 

We have also developed a mean field approximation for the thermodynamic factor in Dc, 
while the mobility is estimated numerically from the Dynamical Mean Field theory. The 
theory reveals that the behavior observed in Dc{0) is due to an entropy- induced repulsive 
interaction. We also extract this interaction numerically from the pair correlation functions 
for athermal chains. It is interesting to note that the diffusive dynamics of polymers in 2D 
is fundamentally different from the 3D case, where entanglement effects dominate in dense 
melts of longer chains |25| . 

The functional dependence of Dc{0) on the coverage has interesting consequences for 
the profiles of spreading films in the submonolayer regime. When the entropy-generated 
repulsive interactions dominate, the droplets assume a spherical cap type of steep shape. 
Contrary to the suggestion of Ref. 0, no assumptions about the film edge being a phase 
boundary between a condensate and vapor need to be evoked here. With increasing at- 
tractive interactions these shapes evolve towars the Gaussian shape. If these interactions 
dominate and Dc(0) is a deacreasing function of 9, profiles emerge that are narrower than 
Gaussians. Thus, the submonolayer spreading experiments constitute a sensitive measure 
on the role of interactions in the diffusive dynamics of polymers on surfaces. 
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FIG. 1. A typical late-time submonolayer spreading profile from experiments of PDMS spread- 
ing on a silver surface || (filled circles). The height of the experimental profile is measured in 
An gstroms (A). Dashed and dotted lines are Gaussian and spherical cap fits, respectively. Solid 
line shows the solution of Eq. (1) using Dc{0) shown in the inset. Dq(0) has been normalized by 
the diffusion coefficient in the zero coverage limit Dq(0). 




FIG. 2. Typical results for the density profile of a droplet from numerical solutions of Eq. (1) 
with different choices of the function Dc{9). 
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FIG. 3. A typical configuration of a chain in the FB model with Nfb = 6. The segments are 
shown by squares. Stars denote the sites blocked by the lowermost segment. An allowed move of 
one of the segments is shown by a dashed line. 
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FIG. 4. Dc{9) for three different values of J as computed from MC simulations of the FB model. 
Circles are for J = 0, squares for J = —0.5, and triangles for J = —1.0. Solid lines are only guides 
to the eye. The curves are normalized by the single monomer diffusion coefficient D\ in the zero 
coverage limit. The error bars are of the size of the symbols, or smaller. 
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FIG. 5. The tracer diffusion coefficient Dt(9) for three different values of J. Circles are for 
J = 0, squares for J = —0.5, and triangles for J = —1.0. Solid lines are only guides to the eye. 
The curves are normalized by the single monomer diffusion coefficient D\. The error bars are 
smaller than the size of the symbols. 
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FIG. 6. A comparison between three experimental submonolayer profiles measured after 50, 80, 
and 150 minutes following deposition of PDMS on silver PJl^] (circles), MC simulations of 2D 
circular droplets from the FB model with Nfb = 6 (squares), and numerical solutions of the 
nonlinear diffusion equation (solid lines). See text for details. 
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FIG. 7. Typical spreading profiles obtained from the MC simulations with (a) J = —0.5 and (b) 
J = —1.0 circles. Solid and dotted lines denote spherical cap and Gaussian fits, respectively. The 
insets show simulated spreading profiles at three different times. 




FIG. 8. (a) Coverage dependence of w\{9) (three upper curves) and W2{9) (three lower curves) 
for J = (circles), J = —0.5 (squares), and J = —1.0 (triangles). Results of fitting [20 are 
shown by solid lines. In the zero coverage limit, both quantities have been computed numerically 
for a single fluctuating chain and thus do not depend on the value of J. (b) Calculated coverage 
dependence of Xq for J = (circles), J = —0.5 (squares), and J = —1.0 (triangles). Solid lines 
show results of numerical calculations using Eqns. (8) and (10). 
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FIG. 9. (a) Coverage dependence of m{9) for J = (circles), J = —0.5 (squares), and J = — 1.0 
(triangles). The corresponding results from the DMF theory are shown by solid, dotted, and 
dashed lines, respectively. The data have been normalized by single segment mobility m\ in the 
zero coverage limit, (b) Comparison between m{9) and Dt{9) with m{9) plotted as in (a) and 
Dt(0) shown for J = (solid line), J = —0.5 (dotted line), and J = —1.0 (dashed line). 
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FIG. 10. Solid line shows the effective potential I4(r)/^T between the athermal chains in 
the FB model (see text for details). The dotted line denotes a typical Lennard- Jones poten- 
tial for comparison. The inset shows seven scaled correlation functions g(f) for Nfb = 6 and 
9 = 0.25, 0.40, 0.50 and 0.75, and for Nfb = 12 and 9 = 0.25, and 0.50, and for N = 48 the case 
9 = 0.25. See text for details. 
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